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Abstract 

We extend our earlier model of qq mesons using relativistic quasipotential (QP) 
wave equations to include open-flavor states and running quark-gluon coupling effects. 
Global fits to meson spectra are achieved with rms deviations from experiment of 
43-50 MeV. We examine in-medium effects through their influence on the confining 
interaction and predict the confining strength at which the masses of certain mesons 
fall below the threshold of their dominant decay channel. 

PACS numbers: ll.lO.Qr, ll.lO.St, 12.40. Qq 
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Interest in calculating the properties of hadronic matter has grown in recent years be- 
cause it is believed that understanding the behavior of these properties in high density and 
temperature environments is necessary for interpreting the results of relativistic heavy ion 
collisions. It is expected that observables such as masses, widths, and couplings will change 
in hot, dense environments like those achieved in experiments attempting to create conditions 
for the phase transition to the quark-gluon plasma [1]. 

Several papers have treated constituent quark models using two-body relativistic wave 
equation treatments [2,3]. We here improve and extend the model of Ref. [3] by including 
open-flavor mesons in the spectrum fits and incorporating a running quark-gluon coupling. 
We then provide initial results for in-medium effects through their influence on the confining 
interaction. 

The interaction kernel consists of a one-gluon exchange interaction in the ladder approx- 
imation, VoGE, and a phenomenological, long-range confinement potential, Vcon- We treat 
the potentials in momentum-space according to the methods described in Ref. [4]. This 
interaction takes the form 



the string tension a is the strength of the confining part of the interaction. We adopt a 
scalar Lorentz structure Vcon since this choice is supported by lattice results [5] and by 
phenomenology, e.g. [3]. 

There are several choices of relativistic two-body wave equations which can be used to 
treat the interactions of Eq. (1). In this investigation, we employ two different equations, 
a spinor version of the Thompson equation [6] and the equation introduced in Ref. [3] . For 
convenience, we will refer to these equations as "reduction A" and "reduction B" respectively. 
These two integral equations result from different choices of three-dimensional propagators 
used to reduce the Bcthc-Salpctcr equation from four to three dimensions. 

Assuming a bound state of quarks with equal masses in the center-of-mass 0~ channel, 
the reduction A equation reads 




Here, ag is the strong coupling, which is weighted by the meson color factor of |, and 



2 



(2a; - E) = ^^h{q, q') + ^^I,{q, q'), (2) 



and the reduction B equation reads 



In each case 



and 



{Aoo' - E^) = ^h{q, q') + —Hq, q')- (3) 

oTiq nq 



h{q,q')= WQo(^) , ^l{q'\ (4) 



hiq, q') = hm|^ H W|Qo(^0 , + \z'Q.{Z') - ll ^ (g')- (5) 

The quantities g and q' are magnitudes of relative three-momenta, 0;*^'-' = [m^ + g*^')^]-^/^, m is 
the constituent quark mass, Qq is the Legendre function of the second kind with argument 
Z — {q^ -\- q''^)/2qq' or Z' — {q^ + g'^ + //^)/2gg', and £^ is the energy to be solved for along 
with the amplitudes ■0° . The '+' subscript on the amplitudes indicates that only the positive 
energy components have been retained in these equations. 

We have improved the calculated spectra of Ref. [3] by incorporating a running strong 
coupling into the model. Specifically, we have constrained the coupling to run as in the 
leading log expression for as, 

where f3i = 11 — 2n//3 and ^/ is the number of quark flavors. The most precise measurements 
of as have been at Mz, the mass of the Z-boson [7]: q;s(a*^ = -^I) — 0.12, and the running 
couphngs used in the following calculations are anchored to this point. We have elected to 
relate to the meson mass scale through 

Q'^eMLson + P', (7) 

where the parameters ^ and p are determined by the fit. The second term, if it is large 
enough, ensures a finite saturation value of as as goes to zero. 
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We have also expanded the work of Ref. [3] by including open- flavor mesons in our fit. 

The constituent masses, the string tension o", and the parameters ^ and p are adjusted 
to minimize the RMS deviation between the calculated meson masses and the experimental 
masses hsted in the Particle Data Tables [8] . 

There is the well debated issue of whether the low-mass pion (and later, when open-flavor 
states are considered, the K meson) should be included in a constituent quark model picture, 
even a relativistic one. We adopt the philosophy that if a model adequately reproduces the 
TT and K masses, then these states may be retained in the calculations using that model. In 
the present case, reduction A does poorly at representing these low-mass mesons, and these 
states are therefore not included in calculations using reduction A. 

The global flts includes 45 states for reduction A and 47 state for reduction B. The 
fltted parameters are given in Table 1. We list the complete calculated spectrum from both 
reductions in Table 2. 

We turn now to a consideration of how this meson model may be used to treat in-medium 
effects on mesons. In-medium effects may be expected to modify primarily the long-range 
part of the interactions. We therefore focus our attention on the conflning potential in 
this simple demonstration of the utility of our meson model. We use reduction B and the 
parameters from the global flt of Table 1 for these examples. 

In order to make our calculations independent of a particular model for the medium 
effects on the string tension, and to avoid the necessity of invoking an equilibrium picture of 
the medium, we plot the meson masses against a/ao, where ctq is the vacuum value of the 
string tension. Given a model for the behavior of string tension with respect to a particular 
variable such as the temperature T, this plot of mass vs. (t/(Tq can be converted into a plot 
of mass vs. temperature. 

It is especially interesting to predict whether meson masses fall below thresholds resulting 
in the closing of decay channels and leading to signiflcant modiflcations in meson widths. 
For this purpose we evaluate the masses of ip" and the D meson as a function of decreasing 
string tension and display the results in Figure 1. 

The calculated -0" and D masses each differ from the measured masses by about 30 MeV. 



4 



In order to make the curves of Figure 1 pass exactly through the and D (x2) measured 
masses for g/gq — 1, we plot M{G)Mexp/ M{gq) on the vertical axis, where M{g) is the 
meson mass calculated using the in-medium string tension cr, and M{gq) is the meson mass 
calculated using the vacuum value of the string tension ctq. 

Figure 1 shows that the ip" mass crosses the DD threshold at a mass of about 3.71 GeV 
when the string tenion falls to about 0.9(To. For string tensions less than this, the avenue for 
■0(3770) decay to DD is blocked. Since this is by far the dominant decay channel for the -0" 
in the vacuum [8], the model predicts significant modifications in its decay signatures will 
result when medium effects decrease the effective string tension by more than about 10%. 
A similar calculation using reduction A also shows this "cross over" at u = O.Qctq. 

In Figure 2 we show similar calculations for the masses of the and the KK threshold 
as a function oi g/gq. The figure shows that the mass crosses the KK threshold at a mass 
of about 9.8 GeV when g/gq ~ 0.67. 

The results displayed in Figs. 1 and 2 may be interpreted within the framework of 
particular models of medium effects on the string tension. For example, several methods 
have been proposed for modeling the temperature dependence of the string tension [9-12] 
and they exhibit the common feature of predicting a monotonically decreasing string tension 
for increasing temperature. 

According to a model for g{T) taken from Ref. [12], as a particular example. 



where Tc is the critical temperature and T is the temperature of the medium. Applying this 
model to the ip" — DD plot in Figure 1, our calculation predicts that the -0" will cross the 
DD threshold at a temperature of 0.43 Tc. 

We may compare our result of 0.43 Tc with the value of 0.56 Tc obtained in the study 
of Ref. [14] which also used Eq. (8) to model g{T). This difference may be attributed to 
the use of different constituent quark models. Ref. [14] used a non-relativistic Schrodinger 
equation treatment to compute the ip" mass, and a bag model to determine the location of 
the DD threshold by using scaling arguments to relate bag pressure to string tension. Our 




(8) 
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model treats all states within a common relativistic framework. 

Turning now to the — KK system of Fig. 2, if we again use Eq. (8) to model the 
temperature dependence of the string tension, we observe that the KK decay channel of the 
(f) will close at a temperature of about 0.74 Tc. 

We may relate this result to that of a previous study [15] which used a Nambu-Jona- 
Lasinio model of in-medium and K masses. In this previous study, the KK decay channel 
for the (f) was found to close at a temperature of 

T ~ 0.45 mi, (9) 

where mi is the value of the light constituent quark mass in a free hadron. Ref. [15] presents 
sample results for light-quark masses in the range 300 MeV < mi < 600 MeV. Using the 
light constituent quark mass obtained in our spectrum fits, mi ~ 340 MeV (see Table 1), Eq. 
(9) predicts that the KK decay channel will become closed for decay at a temperature of 
about 150 MeV. We see that our result of 0.74 Tc would be comparable provided we assumed 
Tc ~ 200 MeV. 

Our relativistic in-medium calculations extend the work of previous non-relativistic po- 
tential models by providing a consistent relativistic framework for evaluating in-medium ef- 
fects, through the dependence of our model parameters on properties of the medium. Work 
in-progress includes calculations of a large set of vacuum and in-medium observables for all 
meson states. 

This work was supported in part by the U.S. Department of Energy under Grant No. DE- 
FG02-87ER40371, Division of High Energy and Nuclear Physics. JPV acknowledges support 
from the Alexander von Humboldt Foundation. We acknowledge valuable discussions with 
J. W. Qiu, S. P. Klevansky, R. Vogt, J. Hill, F. Wohn, and W. Weise. 
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Table 1. 

Best values for the fitted parameters for reductions A, B. Best values for the seven 
parameters were determined by a fit to 45 (47) meson states for reduction A (B). "RMS" is 
the root-mean-squared deviation of the fitted meson masses from measured meson masses. 
Spectra are given in Table 2. 



nib (GeV) 
me (GeV) 
rus (GeV) 
m„ (GeV) 
a (GeV2) 

P (GeV) 
RMS (MeV) 



A 



4.65 
1.37 
0.397 
0.339 
0.233 
0.616 
0.198 
43 



B 



4.68 
1.39 
0.405 
0.346 
0.211 
0.444 
0.187 
50 
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Table 2. 

Fitted spectra for reductions A and B. The seven parameters used to obtain these spectra 
are given in Table 1. 
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Figure 1. 

Plot of the Tp" mass and the DD threshold as a function of ct/cto, where a is the in- medium 
value for the string tension and ctq is the vacuum value of the string tension. Masses were 
calculated with reduction B using the parameters of Table 1. Since using these parameters 
does not yield precisely the experimental vacuum values of the ip" and D masses, the 
vertical axis plots M{a)Mexpl M^ao) so that the curves pass exactly through the 
experimental masses at (t/uq — 1. 
Figure 2. 

Plot of the (f) mass and the KK threshold as a function of a/ao, where a is the in- medium 
value for the string tension and ctq is the vacuum value of the string tension. Masses were 
calculated with reduction B using the parameters of Table 1. Since using these parameters 
does not yield precisely the experimental vacuum values of the and K masses, the 
vertical axis plots M{a)Mej.p/M{ao) so that the curves pass exactly through the 
experimental values at a/ag — 1. 
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